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Abstract 


Despite some advances in prospective renewable energy systems, there 
are still abundant obstacles against rapid development and promotion for 
academic research and industrial applications, including complex systems 
and governed equations, high solution time, and comprehensive quick 
response. Thus, a comprehensive and extendable machine learning-based 
approach is proposed in order to achieve process simplification, quick 
response, less energy consumption, and more precision. The process and 
sequences of the method were presented in detail, it is applied to a 
renewable energy complex problem, and it is utilizable for other similar 
complex problems. In the present article, a general MATLAB code was 
generated for 4E analysis of a parabolic trough solar collector, including 
energy, exergy, economic, and environmental analysis. Proper input data 
was generated with Minitab, increased to a sufficient volume, and pre- 
processing and scaling were performed. Then, through the presented 
machine learning approach, an accurate model was trained in order to 
analyze conveniently, simple prediction and optimization, prompt reaction, 
and implement sensitivity analysis and other requirements. In the wake of 
applying the presented approach, 1143 times faster response and 
calculation resulted as one of the achievements, which made it substantially 
more convenient for the analysis process and optimization. An optimization 


was implemented in a simple way by the trained model by Genetic 
algorithm and revealed the optimum categorical and continuous variables 
during very low calculation, very quick time, and extremely accurate. 
Additionally, a sensitivity analysis and a comprehensive parametric study 
were executed by the trained model in a substantially convenient way for all 
input variables for five objective functions, including energy, exergy, heat 
cost, energy based emission cost, and exergy based emission cost. 


1. Introduction 


Renewable energy is now widely regarded as the most effective way to 
reduce combat global warming, fossil fuel consumption, and survive in a 
more sustainable environment. Renewable energy's quota is 2351 GW or 
thirty-three percent of the total required energy. Based on data that 
reported from the Agency of International Renewable Energy, the capacity 
of renewable increase is still being driven chiefly by wind solar energy 
projects. Solar energy, specially, expanded by 24%, 171GW, in 2018, 
demonstrating the widespread acceptance of its _ utilization in 
various applications [1]. 


One of the important, well-known and promising renewable energy 
technologies is the utilization of parabolic trough solar collectors (PTSC) to 
generate steam for electricity production with an Organic cycle, Rankine 
cycle, or Brayton cycle [2] or to transform solar irradiation into usable 
process heat [3]. The PTSC is made up of three primary components: a 
reflector, an absorber, and a working fluid (WF). The trough includes an 
aluminum surface or polished envelope with a reflectance of 95 or 88.5%, 
respectively [4]. To raise the temperature of the HTF, radiation that strikes 
the surface of mirror is reflected onto an absorber pip that is positioned at 
the focal axis of the parabola. Depending on the application, the HTF, which 
may be water, thermal fluids, molten salt, therminols, or gas pumped 
through the absorber tube, absorbs the energy and exergy from the receiver 
and reaches temperatures of up to 390 °C [4]. 


There are several different PTCS setups that are used for power plants or 
heating and cooling systems. These configurations vary depending on some 
of the impacts, such as the kind of integration, the absorber's composition 
and form, the working fluid, the thermal storage system, and the reflector 
and tracking system types [5,6]. 


The performance and efficiencies have a substantial engagement to the 
parameters as mentioned above. It is confirmed that to approach a high 
performance for PTSC, optimization methods are necessary procedures [6]. 


However, optimizing the various real-time input variable manipulation to 
enhance the efficiencies of PTC systems is one of the primary issues. The 
operating conditions, materials' properties, receiver and HTF type, and 
concentrator and apparatus geometry, among other things, having an 


impact on performance of PTC [7]. Numerous research have centered on 
determining PTC systems' thermodynamic behavior utilizing computational 
models and simulation technologies like machine learning, particularly 
artificial neural networks (ANN) [8]. The experimental systems are limited 
by the high costs of operational and maintenance and are time- 
consuming. The benefits of simplicity, speed, and the capacity to handle 
non-linear and complex interactions between inputs and variables make 
ANN simulation tools and its derivatives sucha promising solution for 
complex data processing [9]. 


Evolutionary Methods (EMs) like the Vector Evaluated Genetic Algorithm 
(VEGA) and Particle Swarm Optimization are two common AI methods for 
resolving Multi Objective Optimization Problems (MOPs) (PSO). After 
intelligent tools are available, the Particle Swarm Optimization (PSO) 
domain gives rise to Multi Objective Swarm Particle Optimization (MOSPO) 
[10]. There are two different forms of multi-objective optimization problems 
(MOPs): linear and nonlinear, as well as convex and non-convex. When all 
objective functions and constraints are linear, linear MOP is taken into 
account, whereas non-linear MOP is taken into account when any objective 
function or constraint function is nonlinear. Convex and feasible regions are 
taken into account while considering convex MOP, and non-convex and 
feasible regions are taken into account when considering non-convex MOP 
[11]. 


A numerical simulation was run to examine the effect of twisted tape 
inserts on the parabolic trough collector's absorber tube's circumferential 
temperature differential reduction and augmentation of heat transmission. 
The maximum circumferential temperature differential will drop by around 
10-45%, the Nusselt number will increase by 1-1.7 times, the friction factor 
will rise by 1-4 times, and this will greatly lessen thermal stresses in the 
absorber tube [12]. 


According to the outcomes of an exergy-economic analysis of a 
commercial parabolic trough collector (PTSC), the exergy efficiency is 
14.1% for the greatest total solar radiation of 665 W/m2 in Tehran, Iran, 
and the exergy dissipation rate is 0.1856 kW/$ [13]. 


By using an economic lens, the impact of climatic variations on a PTSC's 
energy, exergy, and environmental impacts was examined. The PTSC was 
tested in 5 Iranian cities to simulate different climates: Rasht, Shiraz, 
Tehran, Abadan, and Sanandaj. According to the findings of this study, 
Shiraz, which has a Mediterranean climate, has the parabolic solar 
collectors with the best thermal energy efficiency among these cities, up to 
71.97% [14]. 


Some papers tried to assess neural network concepts alongside other analysis 
in the renewable energy field, but the assessments are very limited and they had a 
very brief glance at machine learning next to other analysis. On the other hand, 
high solution, optimization, analysis, and response times are the other obstacles to 


academic research and industrial application. In the current study, a detailed 
machine learning approach is presented comprehensively to overcome the 
hindrances and utilize the procedure for the research and industrial application. 
The sequences of the method are illustrated and described, which it is extendable 
and could be used for other complex renewable energy problems. Herein, the 
analysis and optimization of energy, exergy, heat cost, energy based emission cost, 
and exergy based emission cost are implemented for a complex renewable energy 
system based on the proposed machine learning approach with detailed required 
tools and techniques, which is utilizable for other systems. Table 1 presents the 
study gaps and comparisons with the literature review. 


Table 1: Gap study and comparison of the present article with literature review 
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2. System description 


The schematic diagram, system description, and components of the PTSC are 
illustrated in Figure 1. The components of a related operating system may include 
a PTSC, a pump, a storage tank or heat exchanger, a flow meter, a DC-motor, 
piping, a heat transfer fluid, temperature and pressure sensors, a data logger, and 
a computer. 


The HTF gains thermal heat through an absorber pipe, then the _ high- 
temperature HTF flows to a storage tank to save energy. The energy reservoir 
could be consumed at the same time or whenever needed. The storage tank can 
also derive another HTF through a heat exchanger to other aims, such as providing 
required heat energy for a Rankin Cycle, absorption refrigerant cycle, domestic hot 


water, etc. A computer can detect the thermal properties and hydrodynamic 
parameters of the HTF and environmental data by data logger and temperature 
and pressure sensors. The computer can also control the HTF mass flow rate and 
the angle of the trough to increase efficiencies. Despite A number of parallel and 
series PTSCs, with various range of length and HTFs, could be utilized, the main 
cause is applying a proposed machine learning approach to simplification of the 
complex heat transfer part. Thus, the main concentration is on absorber part of the 
PTSC, which investigated in an abundant detail in following sections. 
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Figure 1: The schematic diagram, system description, and components of the PTSC 


3. Methodology and mathematical analysis 


3.1.1. Optical modeling 


The number of minutes that separate solar time from standard time is 
[20]: 


Tg T=46 (1) 
Where the 7; and T, are the solar time and local time, respectively, L,, is 


the local standard time meridian and L,,, is the site longitude, and E is the 
time equation which is computed by the Following equation [20]: 


i, F=229.2(0.000075+0.001868cos B—0.032077sin B—0.01461544 = (2) 
360(n-1) 
365 


where n is the corresponding day of the year (1 sn S365), and on first of 
January, n =1. 


Where B= 


For a horizontal surface, the angle at sunset is represented as [20]: 


w=cos ' (—tan gtan 6) (3) 


Where g and 6 are, respectively, the deflection and latitude angle. The solar 
deflection angle in degrees could be calculated usingé6 [21]: 


360(284+n) 
365 


ro Pl ee Pressure drop in pip 


6=23.45sin (4) 


The pressure drop of flow through the absorber could be calculated to: 
L V,? 


AP"! paar D Og 


(5) 


With 
4in 
puD, 


v= 2 (6) 


Where AP is pressure drop, L is absorber length, D, is absorber diameter, g 
is gravity constant, V, is mean velocity in absorber tube, m is mass flow and 
p is fluid density. 


The fparey for laminar flow is equal to: 
farey= RR (7) 


Where ¥p2 is Reynolds number of absorber tube. 


And f=fpary/4 for turbulent can be estimated by following equation [22]: 


/D. 
Ss Di6g é 24 2.51 
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Where f is friction factor and € is equivalent roughness which for drawn 
pipe is usually equal to 1.5e-6 m. 


(8) 


3.2. Thermal modeling 


The HCE performance model incorporates the equations and correlations 
required to forecast the terms in the energy balance that associate the 
collector type, HCE circumstances, optical qualities, and environmental 
conditions. It applies an energy balance between the HIF and the 
surroundings. 


By preserving energy balance at each surface of the HCE cross-section, 
the four energy balance equations are established: 


Oi deny = c( Lite Tntot) (9) 
Dy conv =; cond (10) 
ascites = au conr® Taavad* Qos cond * Geena , bracket (11) 
Gu conv? Ga rad Cas cond ( if 2) 
Gu cond + as Soldos — Ge conv + Gye (1 3) 


Where Tite: is outlet temperature of HTF and T;,,.. is inlet temperature of 
HTF. The absorber cross section is showed in Figure 2, where the heat 
fluxes are defined and specified in each surface and related components are 
defined. 
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Glass envelope 
Q12conv: Convection heat transfer from absorber inner surface to the HTF (2-41) 
QssolAbs: SOlar irradiation absorption from the incident solar irradiation to the outer glass envelope surface (s—>5) 


Q3solAbs: SOlar irradiation absorption from the incident solar irradiation to the outer absorber pipe surface (s—»3) 
Q23cond: Conduction heat transfer from the absorber outer surface to the absorber inner surface (3-»2) 


42conv: Convection heat transfer from the absorber outer surface to the envelope inner surface (3-4) 
42rad: Radiation heat transfer from the absorber outer surface to the envelope inner surface (3-44) 


Q45cond: Conduction heat transfer from the envelope inner surface to the envelope outer surface (45) 
Qséconv: Convection heat transfer from the envelope outer surface to the ambient (56) 


Qs7rad: Radiation heat transfer from the envelope outer surface to the sky (5-7) 


Figure 2: Absorber cross section, various zones and temperatures 


3.2.1. Convection Heat Transfer between fluid and 
receiver 


In the balance equations, 1 con is the heat transmission by convection from 
the interior of the absorber pipe to the fluid, and it is also equals to [23]: 


Gs c= Dim =] (14) 
With 

k 
hy=NUpe (15) 


Where D, is the inner diameter of the absorber pipe at point 2, T, is the 
inside surface temperature of the absorber pipe at point 1, T, is the mean 


bulk temperature of the fluid at point 1, Nup, is the Nusselt number based 
on D, at point 1, and k, is the thermal conductance of the fluid at T, of point 
1. 


352.1,1. Modeling of turbulent and transitional 
flow cases 
The following Nusselt number correlation [24] is used to represent the 


convective heat transfer from the absorber to the fluid for turbulent and 
transitional scenarios, Reynolds number greater than 2300: 


£,/8|Rp,—-1000|Pr, \{ Pr,\?"* 
wpa 1+12.7,Vf,/8| Pr??—1 |\ Pr, ae 
with 
f,=(1.8210gjo (8 p»\— 1.64) * (17) 
en Laminar flow cases 


When the flow regime is laminar, Reynolds number less than 2300, the 
Nusselt number will be constant and equal to 4.36 [25]. 


Simca Conduction heat transfer through receiver 
thickness 


The conduction heat transmission through the absorber wall is presented as 
occurring through a hollow cylinder [25]: 


Ce aed er T,- T,)/In |D,/D, | (18) 


Where k,, is the thermal conductance of the absorber at its average 
temperature, T, and T; are the surface temperatures within and outside of 
the absorber, respectively, and D, and D, are the diameters of the absorber. 


Seated: Heat transfer from the absorber to the glass 
pip 

The annulus pressure affects the convection heat transfer process. 
Molecular conduction is the primary method of heat transport at low 
pressures, typically less than 1 torr. Free convection is the mechanism at 
pressures above about 1 torr. The temperature difference between the inner 
glass envelope surface and the inner absorber surface causes the radiation 
heat transfer to occur. By assuming that the glass envelope is opaque to 
infrared radiation and that the surfaces are gray, the calculation for 
radiation heat transfer is made simpler. 


3.2.3.1. Convection Heat Transfer 
3.2.3.1.1. Vacuum in Annulus 


Free-molecular convection is used to transport convection heat between 
the absorber and glass pip when the HCE annulus is under a vacuum of less 
than or around 1 torr [26]: 


4 conv= 1D3.By4|T-T,| (19) 
With 
k 
h = std 
*""[D,/2In [D,/D,|+ ba|D,/ D,+1)] i) 
s\2=a)(9y=5) 
ne 2 a(p+1) 21) 
2.331 B(-20)| T.,4+273.15 
=2:381 B20) T3273.15 = 
P,6| 
Siti dv lhie. Pressure in Annulus 


Natural convection is used to transport heat from the absorber to the glass 
pip even if the HCE annulus is no longer in a vacuum condition of larger 
than 1 torr or so. The following equation is utilized as a result [23]. 

1/4 


2.425 k3,|T3—T,|[Pr Rap;/|0.861+ Pry,|| 


qd conv (23) 
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T,—T,|D: 
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For an ideal gas 
oa 
B= (25) 


Where £B is volumetric thermal expansion coefficient and T;3, is average 
temperature of 7; and T,. 


S3.2cds2- Radiation heat transfer between 
absorber and glass pip 


The following equation is used to determine the radiation heat transfer 
between the absorber and glass pip [25]. 


omD,| 7;-T4| 


Feared“ T1765 4[1—&,)D,/&,D,)] (26) 


-— W 
Where ga is Stefan-Boltzmann constant {5.672 x10°° , €, is emissivity of 
m? K* 


selected Absorber coating, €, is emissivity of glass envelope. 


3.2.4. Conduction heat transfer through the glass 
pip 
Similar to the conduction through the absorber wall discussed in Section 
3.2.2, the conduction heat transfer through the glass pip is calculated using 
the same equation. The temperature distribution is taken to be linear, much 
like in the absorber example. With a value of 1.04, the thermal conductance 
is considered to be constant [23]. 


Dede: Heat transfer from the glass pip to the 
ambient 


Convection and radiation are the two way used to transfer heat between the 
environment and the glass enclosure. The convection will be forced or 
natural based on whether there is wind or not. Heat is lost by radiation 
when there is a temperature differential between the sky and the glass 
enclosure. 


3.2.5.1. Convection heat transfer 


A primary source of heat loss is the convection heat transfer from the glass 
envelope to the atmosphere, particularly when there is a breeze. 


Gsene= Neg D| T=, (27) 
k 
hsge=-N Ups (28) 


5 

3.2.5.1.1. No wind condition 
Natural convection will be used to transport convection heat from the glass 
pip to the environment if there is no wind or just very little wind [25]: 
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Prsg=— (32) 
56 
Where Pr;, is the Prandtl number for air at T3,, V5, is the kinematic viscosity 
for air at Ts,, Tsgis the mean film temperature at (7;+T7,)/2, and Rap; is the 
Rayleigh number for air based on the outer diameter of the glass envelope, 
the gravitational constant, and the thermal diffusivity for air at T;,. This 


association holds true for 10°<Ra,,<10"* values. 
3.2.5.1.2. Wind condition 


Forced convection will be used to transfer heat from the glass envelope to 
the environment if there is wind with a wind speed of more than 0.1 m/s or 
such [25]. 


1/4 
Pr, 


Pr. 


For Pr=10 n = 0.36, and for Prs10 n = 0.37. This link holds true for the 
values of 0.7<Pr,<500 and 1<Re,<10°. Except for Pr, which is assessed at 
the glass pip outer surface temperature, all fluid parameters are estimated 
at the ambient temperature, 73. 


Nup;=CRp5Pr¢ (33) 


3.2.5.2. Radiation heat transfer between glass pip and 
ambient 


The temperature differential between the glass enclosure and the sky is 
what causes the radiation transmission [25]. 


O57 rea= 9 D, 1é,| T3— T5| (34) 
Where o is Stefan-Boltzmann constant, €, is emissivity of the outer surface 


of glass pip and T; is effective sky temperature. It is estimated that the 
effective sky temperature is 8°C below the surrounding temperature [27]. 


3.2.6. Solar irradiation absorption 
3.2.6.1. Solar irradiation absorption in the glass pip 
Is sotabs = Wei Nenv Cenv (35) 
with 
Nenv= 1 £2 €3 E4 &5 & Park (36) 


Where q,, is solar radiation per receiver length, 1).,y is efficiency of effective 
optical at the glass envelope and @,,,, is the absorptance glass envelope. E, is 
HCE Shadowing, «, is Tracking Error, ¢,, is Geometry Error, €, is Dirt on 


Mirrors, e is Dirt on HCE, é. is Unaccounted errors, P,; is Clean Mirror 
Reflectance that are equal to 0.974, 0.994, 0.98, reflectivity /p,, , [1+ €,)/2, 
0.96, 0.935, respectively [27]. K is incident angle modifier which is equal to: 


K=cos (6)+0.000884 6—0.00005369 &” (37) 
Where 6 is tilted angle of solar collector. 
3.2.6.2. Solar Irradiation Absorption in the Absorber 


as SolAbs — ei Mianet (38) 
With 
TNi= Nene Toor’ (39) 


Where 77, is the effective optical efficiency at the glass pip, q,, is the 
absorber's absorptance, Tepy is the transmittance of the glass envelope, and 
Neny iS the absorber's absorptance. 


ie ea a Heat Loss through Support Bracket 


Support brackets that extend from the collector structure to the receiver 
support the HCEs at the collector focal line. Every HCE has a support 
bracket at the ends (every 4 meters of receiver length), and the current 
work is supported by two of these brackets. By modeling the support 
bracket as an infinite fin with a base temperature 10 degrees lower than the 
outer absorber surface temperature T3 at the attachment point, the bracket 
losses may be roughly calculated. [25]. 


ne Vh,P,k, Acs, | Tes De Les (40) 


Where, in accordance with section 3.2.5.1, A, is the average convection 
coefficient of the bracket that is dependent on wind speed. P, is the 
bracket's perimeter, k, is its conduction coefficient, A,,, is its minimum 
cross-sectional area, T,,,. is the temperature at the bracket's base, 7; is the 
surrounding air temperature, and Lyc, is the length of the HCE. The HCE 
support bracket has a circumference of 0.2032 m, a cross-sectional area of 
1.613e-4 m?, and a conduction coefficient for plain carbon steel that is 
constant at 48 W/m-K. 


3.3. Exergy analysis 
The outcome of applying the second law of thermodynamic on PTSC is: 


EP Nvenpeene, omit 4 (41) 


Where Fx; .,, is exergy input to the fluid flow, Ex, is exergy flow on entering 
solar radiation exergy flow on incoming solar radiation, Ex, out 1S Output 
exergy of the fluid flow, Fx, is the exergy loss and EX tes is the exergy 


destruction. Quantity of energy flow based on received solar radiation can 
be estimated by Petela-Landsberg-Press equation [28]: 
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The Petela-Landsberg-Press equation is constraint to fully concentrated 
solar radiation, 46200 or so, and blackbody absorbers. Since the PTSC have 
smaller concentration ratio and their absorber is not blackbody,so prior 
Equation is not valid on this condition. Eventually, following equation is 
proposed to estimate the amount of exergy flow that is incoming in PTSC 
from the sun [29]. 
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Where T,,,, is the temperature of sun surface and supposed by 5780 (k) and 
Jis the interaction factor and calculated by following Equation [30]. 


0.25 
fy, XE_Xa 


l= 
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(44) 


Where fy, is the view factor that takes the geometrical relation between the 
radiation source of the absorber of radiation, €, is the dilution factor that 
hits on a surface from a radiation source of view factor and f, is Geometric 
Factor of the emitted radiation by the Absorber. 


Calculating the usable exergy involves comparing the exergy of the fluid 
imported into and exported from the collector [31]: 


AP 


LEx,=Ex, .,-Ex AX 
. eine £,€6=O,—1,X C,X Tomy X10 Fat nx Tome p,X T; (45) 
Where AP is the pressure drop value. 
The value of exergy loss is equal to [32]: 
Ex = Ex, optt Ex, thermal (46) 


Where FE X) op: iS the optical exergy loss, and E X) tnerma) 1S the thermal exergy 
loss, which are calculated by two previous Equations, respectively [33]. 


FG p= la | REX (47) 
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HB 


1-7 (48) 


The irreversibility that is developed in the heat transfer term is revealed by 
exergy destruction. This parameter specifically denotes a potential work 
that vanishes when thermal energy is moved from a hot to a cold resource. 
The extent of energy destruction is evident [34]: 


BK j= EX geese? BE Mijes oy (49) 


In PTSCs, the exergy destruction contains two main terms; one is between 
the Sun and the absorber tube, and the other is between the absorber and 
heat transfer fluid [32]. 
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Eventually, the exergy efficiency can be calculated from following Equation 
[34]: 


_ Ex, (52) 
Nexergy ae Ex, 


3.4. Economic modeling 


The total net present cost (NPC) of PTSC is given by following equation 
[35k 


NPC= CRF x(Cii OA M+ C,)é (53) 
With 
CRF=ii, (54) 


Where CRF is capital recovery factor, Co,y is operation and maintenance 
costs which is equivalent to 2% of the initial investment cost, C, ($) is the 
cost purchase, i is interest rate which is equal to 2% and n is number of 
years. 


The heat production cost of PTSC is given by following equation [35]: 
_ NPCIN 
Tu 


Where N is annual number of operation hours and q,, is useful heat transfer 
to the HTF. 


The purchase costs of the proposed system components are presented in 
Table 2 [35,36]. 


Cy (55) 


Table 2: The estimated system components' purchasing prices 


Numb refere 
~ Components Unit Price aee 
1 Solar collector $im* 148 [37] 
Heat Transfer fluid 2 
2 (water) sim 4 [37] 
3 Storage tank $lm? 27 [37] 
A Pump component $ 3540|/ We)" [36] 
O.0% Environmental analysis 
3.5.1. Energoenvironmental 


Environment analysis of a system was introduced by Faizal et al. [38] 
according to following equation: 


Xco2 = Vco2 x Q, x Psat (56) 


Where X,,9 is the released CO, amount at the considered time in kgCO,, and Y,o2 
is the CO, emission value associated to the reference energy system gained by the 
life-cycle assessment approach. Additionally, Q, is the generated power rate of the 
reference system, and ly,o:xing iS the system's working time. 


The life cycle assessment, LCA, is one of the environmental assessment techniques 
that is now used as a common and all-purpose technique in the evaluation of the 
environmental impact of processes and goods. The life cycle assessment is the 
criterion for determining how each action will affect the environment. The third 
component of a sustainable assessment is completed by this parameter together 
with economic and technical evaluations. Along with the economic and 
technological effects, it is important to consider the environmental ones. This 
technique evaluates the environmental effects of each human activity or product, 
including the processes involved in its extraction, manufacture, lifetime, disposal, 
and recycling. Table 3 presents the comparison among different values related to 
carbon dioxide emission value in various systems [39]. It should be to state that the 
values showed in Table 3 are obtained by the LCA method. 


Table 3: Different values related to carbon dioxide emission value in various systems 


carbon dioxide 
emission amount 
[41] (kg CO2/kWh) 


carbon dioxide emission 


system type = amount [40] (kg CO2/kWh) 


Hydraulic 0.0037-0.237 0.010-0.013 
Wind 0.0097-0.1237 0.009-0.010 
Solar thermal 0.0136-0.202 0.013 
Nuclear 0.0242 0.066 
Biomass 0.035-0.178 0.015-0.041 
Solar PV 0.0534-0.250 0.032 
Gas 0.6076 - 
Oil 0.7421 0:778 
Coal 0.9753 0.960-1.050 
Biogas - 0.011 
Natural gas - 0.443 
Geothermal - 0.038 
Fuel cell - 0.664 
Solar collector 
[42] 0.00647 
Dedede Exergoenvironmental analysis 


Another approach that suggested for a system’s environmental assessing is 
exergoenvironmental (EXEN) analysis. The analysis is shown as following 
equation by Caliskan et al. [42]: 


X, 


ex,co2 


A system's environmental analysis is comparable to its EXEN analysis. The 
difference is that energy values are used in environmental analysis systems, 
whereas exergy values are used in EXEN analysis. Because the values 
derived from exergoenvironmental analysis are lower than those derived 
from energoenvironmental analysis, the energoenvironmental analysis is, in 
general, a more wary criteria. 


= eth eg ene (57) 


3.6. Enviroeconomic analysis 


3.6.1. Energoenviroeconomic analysis 


According to following equation the enviroenconomic analysis of solar 
systems has been presented by Deniz et al. [43]: 


Co02=Xco2 ¥ Cco2 (58) 


Where C,,. is the cost of carbon dioxide emissions and C,,,. is the environmental 
economic parameter. The cost of CO, emissions generated in the pertinent system 
is, in fact, the parameter that was obtained. It can also be viewed as a technique 
for assessing how energy systems affect the environment [44]. Several techniques 
and strategies can be used to reduce CO, emissions as much as feasible in order to 
prevent climatic variation similar to some events, such as global warming. The 
carbon released cost technique with an economic perspective aims to increase the 
relatability of the environmental effects. This technique combines punitive and 


incentive strategies to lessen the negative effects of energy systems on the 
environment. 


Based on varied legislation, the cost of the generated carbon varies by country and 
ranges from 13 to 16 dollars per tonne of carbon dioxide [45]. This analysis takes 
into account a mean value for price changes of 14.5 ($/ton), which is equal to 
0.0145 ($/kg) of carbon dioxide [46]. 


3.6.2. ExergoEnviroEconomic analysis 


An analysis of a system's energy and economics is similar to its analysis of its 
environment and economy. The calculations in the exergoeconomic analysis are 
based on the exergy parameter, whereas those in the enviroeconomic analysis are 
based on the energy parameter. In this regard, the exergoeconomic analysis for 
solar energy systems is suggested using Equation [46]. 


The interactions between the six broad parameters that have been introduced are 
shown in Figure 3, together with the contributions made to each of the four 
quantities energy, exergy, environmental impact, and economic domain. The 
connection between the exergy and environmental domains widens the EXEN area, 
as seen in Figure 3. The ENEN domain develops as a result of interactions between 
the energy and environmental domains. EXENEX is produced in a similar manner 
through the intersection of the energy, environmental, and economic realms. The 
interplay between the energy, environmental, and economic realms leads to the 
development of ENENEC. 


Environmental Analysis 


ES > Xozimp (2) 


Enviroeconomic /{ 2&0 
Analysis BEES 


Money Policies 


Figure 3: Detailed interaction between environment and economy 


eZ Machine Learning 


A group of methods and tools known as artificial intelligence can be used 
to empower computers to carry out difficult jobs. The first artificial 
intelligence systems were created in 1950, but for a very long time there 
were no practical applications because computers could not process 
massive amounts of data. A breakthrough in artificial intelligence has 


occurred since 2010 [47]. Since there have been so many advances in the 
field of computer science, anyone can now create their own artificial 
intelligence. 


As it can be seen in Figure 4, neural networks are a sub-section of 
machine learning approaches, and Artificial Intelligence and Data Science 
are in common in Machine learning. Machine Learning includes, Neural 
Networks, Evolutionary Algorithms, Search and Reinforcement learning. 
This progress has brought up both new opportunities and risks. Cars can 
crash, patients can be wrongly diagnosed, people's privacy can be invaded, 
and people can be unfairly stigmatized based on their ethnicity and 
appearance as a result of errors made by AI algorithms. Therefore, it is 
crucial to inform anybody working with AI about what AI is, how it 
functions, and what errors can occur while utilizing artificial intelligence. 


Artificial Intelligence 


Data science 


Machine learning 


Neural 


networks Search 
Linear 
Rein ment 
learning regression 


Figure 4: Relatives among, Machine Learning, Artificial Intelligence and Data Science, Neural 
Networks and other related fields [47] 


3.8. Neural network algorithm 


Artificial Neural Networks (ANN) have had a striking expansion in 
recent years, both in basic and applied research. This is due to the fact that 
ANN is a cutting-edge approach that has never been used before and can 
resolve challenging linear and non-linear issues for which traditional 
mathematics, algorithms, and methodologies have been unable to come up 
with an adequate and acceptable answer [48]. 


Machine learning algorithms known as ANNs are based on how the 
human brain functions, which is included several artificial neurons related 
to a parallel network. These networks function by connecting the inputs and 
outputs of the entire system using non-linear mapping techniques [49]. 


The fundamental structure of an ANN consists of an input layer, one 
or more hidden layers, and an output layer. Each neuron in the output layer 


receives the information from the layer above. Each output neuron receives 
information from the buried layer. Each neuron in the hidden layer receives 
information from the input layer, whereas the input layer receives 
information from the experiment. The link between neurons and layers is 
connected with a weight value and a compensation value known as bias, 
allowing them to cooperate. The mathematical equation that describes how 
the network works is as follows: 


n= DenlIW, p Pt Dy, 


(59) 


IW, , are the coefficients of the connection weights between the input 
layer and the hidden layer, P, are the input variables, and 5, is the bias by 
corresponding to the neuron in the hidden layer, where K is the number of 
input data (k=1,2,3---K). The data is stored as a set of biases and weights 
for connections. Following the addition of the weighted inputs from various 
connections and the compensation value, a transfer function is applied to 
neuron J of the hidden layer: 


K 
aj=F| > OIW, ,P.+ Dy. (60) 
k=1 
In which Y, experiences the same process as described by: 
Ss 
Y,=F| >) O(LW, a7 D,,| (61) 
s=1 


Transfer functions are used to return output for a specific range of 
input values in the hidden layer. The most widely used transfer functions 
are the tangential-sigmoidal (TANSIG) and logarithmic-sigmoidal (LOGSIG), 
which are determined by Transfer function using a logarithmic sigmoid: 


F\n))=LOGSIG |n,|= reg O<F last (62) 
Hyperbolic tangent sigmoid transfer function is given by: 
1-e°” 2 
F\n,|=TANSIG | n |= = -1,-1<Fin)<1 63 
| | | | lae Ov leeC% | | ( ) 


To enhance experimental data prediction in the output layer, the 
linear transfer function (PURELIN), which is defined as: 


Linear transfer function: 


Ss 


VAL W, at dy»| (64) 


s=1 


DOLW, ,aj+.By |] ¢ PURELIN 
s=1 


é 
It is important to note that before training the neural network, 
standardizing each input variable is a valuable tool. Interval operating of 


non-linear transfer functions is the first step in training for increased 
efficiency. The input normalization technique typically scales the data using 
a linear mapping known as the min max scale. 


The following equation is utilized to standardize the input variables 


[50]: 
PP) 
Die sp i Vista es Ye (65) 


Where P;y represents the normalized variable, P, represents the 
variable's actual value prior to normalization, P,,,, and P,,;,, express the 
variable's maximum and minimum values prior to normalization, and 
Y max’ \min State the variable's maximum and minimum values following 
normalization. 


To evaluate the effectiveness of the ANN model, experimental data 
were compared to simulated data using statistical analysis. Three often 
utilized test parameters are the Root Mean Square Error (RMSE), Mean 
Absolute Percentage Error (MAPE), and Coefficient of Determination |R?| 
[51]. Ideal models are those with |R?| values that are near to 1, as well as 
MAPE and RMSE values that are close to 0. To verify the linearity and 
accuracy of the model, the intercept-slope test was also carried out with 
slope = 1 and intercept = 0. 


2 
pase) 22% 7 Xout ~ «| (66) 
0 Xout, exp (4) Xouts,sim(i) 
= Xout exp (4) (67) 
MAPE= _ x 100(%) 
= 2 
pD [Xout ,@Xp (i) Xout ~(s| 
R’?=1-— i=1 = (68) 
Y Deed 
1=1 
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where Xjut exe 1 Xoutexp (i) Xout ,~ (i) iS the gained output value by ANN, and 


Xoutexp (i) iS the experimental output value. 


The matrix of weights was used in a sensitivity analysis after the ANN 
model had been trained to determine the relative contribution of each input 
variable to the output variable and to determine bias. According to the split 
of the connection weights, Garson [52] proposed the following equation: 


(69) 


k=N, m=N, tial - 


Where I, represents the relative weight of the j” input variable on the 
output variable, and N; and N, represent the number of neurons in the 
hidden and input layers, respectively. The input, hidden, and output layers 
are indicated by the superscripts J,h, and o. Input, hidden, and output 
neurons are identified by the subscripts k,m, and n, respectively. 


3.9. Proposed Machine Learning Approach 


Solving a striking number of equations in engineering requires high 
time, energy, cost, and bulk of calculation because there are usually 
multiple equations with multiple unknowns, multiple physics, complex 
physics, non-linear equations, or coupled equations. Sometimes, solving 
equations and extracting solutions to respond quickly and satisfy objectives 
is an essential requirement. The approach with these hindrances is not only 
executable with elevated hardware computers, but also difficult for custom 
personal computers. Thus, it is being tried to solve these barriers in a subtly 
affordable way. In this respect, applying a confirmed and accurate machine 
learning model could be troubleshooting and redeeming. Figure 5 illustrates 
the applied approach as a comprehensive depiction, including procedure 
selection, pre-processing of data, post-processing of data, training, desired 
circumstances, and achievements. 


In addition, applying an appropriate data scale could have a substantial 
impact on the precision of the model, stability, performance, and time 
period of training. The examples in a training dataset help the models learn 
how to map inputs to outputs. In order to estimates of error on the training 


dataset, an optimization algorithm updates the machine learning model's 
weights, which are initially initialized to small random values. 


The bulk of the inputs and outputs which utilized to train the model is a 
crucial consideration due to the modest weights used in the applied model 
and the utilization of error between predictions and expected values. 
Unscaled target variables on regression issues could lead to exploding 
gradients that cause the training process to collapsel, but unscaled input 
variables might cause a delayed or unstable learning process. 


Is the problem complex or does 
the main function take a long time 
or is a quick response required? 


Analysis and 
optimization by the 
trained model 


Apply the novel presented 
approach in the present article or 
not? 


Custom analysis and 
optimization 


Generating proper input data appropriate to 
realistic bounds (e.g. by Minitab) 


Undesired circumstance 


Increase input data to a thorough covering 
and sufficient data availability with an 
appropriate programming language (e.g. 


MATLAB or Python) Achieving novel industry and 


research accomplishments, 
investing in new instruments and 
methods, less cost, less energy, 
save time, more accuracy, and 
easy development 


Applying non-linear constraints to excluding the data that 

cause contrast in reality and the fundamental physics of 

the problem (this process is a subset of running the main 
function and through extracting the output data) 


Extracting output data associated with the input data 


through the main function of the proble A facilitator and utilizable 
appliance provides a prompt 
reacting, comprehensive solution, 
and forecasting for programming, 
research, and industry 


Excluding failed output data and applying a 
scaling method in a proper related range for 
input and output, in order to increase stabililty, ’ 
decrease training time, and avoid exploding requirements 

gradients 


An obtained trained model with sufficient 
Training model through utilizing a proper accuracy and quick respondent 


Machine Learning algorithm (the simplified main function) 


Figure 5: Sequences of proposed Machine Learning approach for complex engineering problems to 
simplification and prompt response 


4. Results and discussion 
4.1. Model description, PTSC_ specification and _ climate 
condition 


The illustrated system in Figure 1 have been modeled based the presented 
Machine Learning approach and optimized for Tehran city, capital of Iran, 
which located at longitude +52.5° and latitude +35.7° [53]. 


Table 4 presents the average climate condition in Tehran, such as average 
air temperature, wind velocity, DNI, average monthly hours of sunshine and 
average monthly solar incidence angle. For this system to be applicable 
along the days of the year, average data considered for this optimization. 


Table 6 shows the specifications of PTSC, and The heat transfer constants 
and coefficients for each annulus gas are listed in Table 7, which are 
studied in this article. One main code is generated in MATLAB software to 
calculate the parameters of the PTSC, one code is prepared for the GA 
algorithm, and one code is written to considering non-linear constraints of 
the PTSC optimization. 


4.2. Validation of model 


The modeling and generated codes in this work are validated with Ref [54]. 
Climate circumstances examined in reference [54], as well as the PTSC 
specification, are regarded as input data. Figure 3.2, and Table 3.3 of that 
reference are considered. In those tables and figures outlet temperature of 
HTF, the difference between the average HTF and the ambient temperature 
( Taira), useful heat transfer and energy efficiency are calculated for the 
mean values along the year. Table 5 shows how the outcomes of this 
modeling were compared to reference experimental data [54]. 


Table 4: The average climate condition in Tehran including Air temperature, wind velocity, DNI, 
average Monthly Hours of Sunshine, and average monthly solar incidence optimum angle 


month 


Value 
S 
| average 
Air Wind Average monthly 
temperat velocit DNI ieeney solar 

es ours 0 ee 

’ [KWh/n Sunshine peenate 
om an sti 


Jan 5.8 2.6 Dad Tak 63.1 


Feb 8.4 2.6 3.5 7 53.8 


Mar 16.9 4.2 4.5 9.3 38.4 
Apr 16.7 3.2 5.5 r.2 17.5 
May 21.3 3.6 6.4 11,8 0.1 
Jun 29.1 3 79 12.2 26.1 
Jul 33.9 2.9 7.0 12 -6.1 
Aug 31.6 2.9 6.5 12 7.5 
Sep 26.9 2.9 5.6 11.6 26.5 
Oct 19.1 2.8 4.0 7.5 43.5 
Nov 10.9 Disk Zod ce 54.0 
Dec 8.9 Ly 2.4 vee) 57.1 
Annua 

wvera 1913288 4.BE GS 28.9 
ge 


Table 5: Comparison of presented modeling outcomes with experimental data [54] 


Parameter Reference [54] Present model Error (%) 
Toutiet( °C) 124 123.7 0.2 
T ee 91.9 91.6 0.3 
Heat gain (w/m) 3402 3206 4.8 
Energy efficiency 72.5 71.85 0.8 
(%) 


The maximum error of the results is 4.8 % which is due to determine the 
variation in heat loss between theoretical and experimental states. 


Table 6: The specifications of PTSC [54] 


lead Parameter Unit Value 
1 CO snsorber = 5 0 9 5 5 
2 e absorber _ é 0) ‘ 1 5 
3 k fees W/mK 1.04 
4 E ies —%6 0.86 
5 Tolass —é 0.965 
6 O vas —%é 0.02 
Mirror 
7 reflectively “i 0.935 


Table 7: Heat transfer coefficients and constants for each annulus gas 


Annulus 
Gas 


4.3. Application of proposed Machine Learning 
Approach 
In the present article, the comprehensive proposed procedure is implemented 


for a complex problem that deals with energy, environmental, and economic issues 
and their consequences. This approach is defined and performed with all the 


details and is extendable to the other complex problems in engineering, including 
heat transfer, thermodynamics, fluid mechanics, and any other disciplines of 
science which are collided with the mentioned obstacles or any other complexities. 
Dealing with non-linear multiple equations and multiple unknowns, a few times 
and in different places, is one of the complexities of this article’s problem. Solving 
non-linear five equations with five unknowns for extracting temperature 
distributions in various zones causes to consume high time, energy and amount of 
calculation. For this case study, the period of time needed for each main function is 
equal to 8 seconds, and optimization period time is estimated to: 


250| generations) x200( populations) x8 (runtime of the main function)=111. 


These periods are too long and energy-consuming. Additionally, they do 
not provide acceptable time solutions to theoretical problems, and do not 
satisfy quick responses in industrial applications. Therefore, an attempt to 
simplify the issues and overcome the obstacles is an essential requirement. 


The primary standard input data was generated with Minitab by creating a 
central composite full of response surface design. After that, the number of 
continuous and categorical factors was defined as seven and three, 
respectively. Then, related upper and lower bounds were assigned. After 
creating the adjustment options, 6912 base runs were generated as input 
data. Then, a MATLAB code is generated to increase the input data to 
provide sufficient data. The input data was increased to 173631 data. Next, 
the main function code of the system ran for the number of data. While 
running the main function, the essential point is to apply the non-linear 
constraints to avoid non-physical data. There is some failed output data as a 
result of satisfying the non-linear constraint. After skimming the failed data, 
the bulk of utilizable data remained for 45835. 


The following content includes training a proper and accurate machine 
learning model to achieve immediate response, consuming less time 
process, less calculation requirements, and energy due to treating the 
problem conveniently. In addition, the effect of various parameters of 
training is discussed. 


For training of the model of PTSC, several inputs and outputs are 
considered. Outputs include energy efficiency (%), exergy efficiency (%), 
and heat cost ($/kwh). The heat cost function includes the cost function ($) 
and useful heat gain (W), which is related to energy efficiency. According to 
the changing of seven continuous variable values with related bounds in 
Table 8 and three categorical variables in Table 9, the target is to provide 
an accurate model between inputs and outputs. It should be noted that 
some of the outputs are not independent, and they could be calculated from 
the other independent outputs. It is possible to consider more outputs to 
elevate the accuracy of prediction or based on desired parameters by 
presented tips in methodology section. The investigated inputs and output 
functions are presented in Table 10. In addition, the related non-linear 


1Ar 


constraints are so vital to exact and realistic prediction while training. The 
non-linear constraints are presented in Table 11. 


Table 8: The variation of seven continuous variable values with related bounds 


- Variable a Lower bound UPPEr 
O. m. bound 
1 Dimenansorber 12 0.01 0.08 
2 R,-R, m 0.005 0.05 
3 £ m 2 100 
4 W m 1 6 
5 m kg/s 0.0001 12 
6 Pannulus Da 0 88223 
7 ; °C 20 80 
outer absorber 1 D tecinotes = i ieeaneothen +6 
0.003 
Deutergiass Mm — Dayutergiass = Dinner giass *9-003 


Table 9: The three categorical variables 


Variabl 
e 
Thermin : . : Heat 
Sylther Thermin Thermin Thermin 
ol72  ms00 olTpDi2 Water olVP1 ol66 ‘anster 
fluid 

Aluminum, Absorbe 

2024, r tube 
Temper- ence aie materia 


T351 


Hydroge ie Annulus 


Argon i Gas 
Table 10: Investigated inputs and outputs for training 
Sign Dim. Description 
Design D. a Inner diameter of absorber 
variables ee eeroevon tube 
RR; - Distance between absorber and 
envelope 
m kg/s Mass flow rate 
T; °C Input fluid temperature 
L m Collector length 
W m Collector aperture width 
— pa Pressure between absorber and 
envelope 
XA shsatbar — Absorber tube material 
Aap — Heat transfer fluid type 
». fee — Annulus gas type 
Outputs FArst= Nonergy % energy efficiency 
SeCONA= Nexergy % exergy efficiency 
Third=C,, $/kWh heat production cost 


Non-linear constraints 


Table 11: Applied non-linear constraints for main function for training at the same time 


Acceptable bound 


HTF velocity limit 


Velocity <3 m/s 


HTF temperature Tire TS Tat 
HTF pressure Purp= Peat 


Exergy efficiency Nee e 


Figure 6 depicts the trend of training by training data, validation data, 
test data and the best point for each epoch. As it is shown, Mean Squared 
Error (MSE) or Residual Sum of Squared (RSS) is decreasing, and the best 
validation point is at epoch 97 with performance of 2.7263 x10’. 


Best Validation Performance is 2.7263e-07 at epoch 97 
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Figure 6: The variation of Mean Squared Error by Epochs with Training, validation, test data 


Figure 7 illustrates the Error histogram around zero error line for 
instances with 20 Bins for training data, validation data and test data. The 
error, which is equals to targets minus outputs, is distributed between an 
acceptable area, which is between 0.001717 and —0.0017. 


«104 Error Histogram with 20 Bins 
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Zero Error 


Instances 


01025 


-0.02735 
-0.02564 
-0.02393 
-0.02222 
-0.02051 
-0.0188 
-0.01709 
-0.01538 
0.01367 
0.01196 
-0 
-0.00854 
0.00683 
-0.00512 
-0.00341 
-0.0017 
7.31e-06 
0.001717 


Errors = Targets - Outputs 


0.003427 
0.005137 


Figure 7: The Error histogram around zero error line for instances 


Figure 8 shows the regression results for training data, validation 
data, test data and all of them in a comprehensive depiction. Outcomes 
confirm that an ideal regression is resulted, which the regression number is 
equals to 0.9999 for test data. In addition to utilize various procedures to 
increase accuracy, this precision is obtained from two main reason, 
including scaling the input and output data to an appropriate range and 
applying non-linear constants to ensure being physical meaningful data. The 
scaling data avoids generating large weight values from training model. 
Large weight values are frequently indicative of an unstable model, which 
may exhibit poor learning performance and sensitivity to input values, 
leading to a larger generalization error. 


Assessing the results under scrutiny emerges the lack of data for 
regression for area that non-linear constraints banned producing data. 
Thus, there are no data in this vicinity, because this data is in contrary to 
reality and the fundamental physics of the problem. Following this approach 
is a substantial step to avoid misdirection, inaccurate machine learning 
training and generating unrealistic outputs from trained model. 
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Figure 8: The Regression results for training data, validation data, test data 


4.4. Optimization of PTSC 


In the wake of providing an accurate proper machine learning model, 
executing optimization for the prior complex problem, is simple, quick and 
convenient for various decision variables to optimization and analysis of 
desired objectives. 


For optimization of PTSC with GA algorithm, by trained model, two 
objective functions are evaluated: exergy efficiency (%) and heat cost 


($/kwh). The heat cost function is included cost function ($) and useful heat 
gain (W) which is related to energy efficiency. According to the changing of 
seven continuous variable values with related bounds in Table 12 and three 
categorical variables in Table 13, the goal is to maximize energy efficiency 
and reduce the cost of providing heat. Table 14 presents the variables and 
objective functions that were examined. 


Table 12: Lower bounds and upper bounds of seven continuous variable values for optimization 


oe Variable es Lower bound SDDer 
er m. bound 
1 Donerabsorver I 0.01 0.08 
a R,-R, m 0.005 0.05 
=) i m Z 4 
4 W m 1 3 
5 m kgls 0.0001 12 
6 Proms ‘Pa 0 88223 
7 i °C 25 80 
DO citer pasion. 2 DP sities sitnay =O misnanenee 
0.003 
outerglass MM Dygegincs =P reverting 40.003 
Table 13: categorical variables for optimization 
Variabl 
e 
ol72, Solther Thermin water Thepmin Therm. tranofor 


Aluminum, Copper 321H Absorbe 
2024, r tube 


Temper- materia 
T351 ] 
Argon Hydroge hae Annulus 
n Gas 


Table 14: Determined variables and objective functions for optimization 


Sign Pen Description 
on 
Design ee Inner diameter of absorber 
variables mneubyenien tube 
Distance between absorber 
R,-R, m 
and envelope 
m kg/s Mass flow rate 
ie G Input fluid temperature 
i, m Collector length 
W m Collector aperture width 
Pressure between absorber 
annulus pa 
and envelope 
pete — Absorber tube material 
HTF — Heat transfer fluid type 
empl — Annulus gas type 
Objective Te " = 
fon. ODj.1"=Nexergy % exergy efficiency 


Oj 2"°=C., $lkWh heat production cost 


The main parameters of Genetic algorithm solver are presented in Table 15. 
Furthermore, the prior applied non-linear constraints should be applied 


again for optimization as presented in Table 12. One of the benefits of 
applying proposed machine learning approach is to decrease optimization 
process time which is shorted to a logical and affordable range and 
descended to 1143 times less. 


Table 15: Main defined parameters of Genetic algorithm solver 


Parameters Multi objective Genetic Algorithm 
Decision 9 
variables 
Size of . 200 
population 
Fraction of 0.65 
crossover 
aaa adaptfeasible 
unction 
Sire intermediate 
unction 
Selection 
é tournament 
function 
Function 10-2 
tolerance 
Number of 250 
generation 
Max stall AO 
generations 
Running time 250 x200 x0.007(second)=5.8 mins 


Table 16: Non-linear constraints for optimization 


Non-linear constraints Bounds 


HTF velocity limit Velocity <3 m/s 


HTF temperature Tupr\ 17S Tat 
HTF pressure Purp= Peat 
Exergy efficiency Reteray =O 


Pareto front 


0.01546 


0.01544 


0.01542 


0.0154 
Optimal Point 


Heat cost [$/kwh] 


0.01538 


0.01536 


0.01534 : ; 
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Figure 9: Pareto front of two objective functions derived from the trained model 


Figure 9 illustrates the Pareto front of these two objective functions 
derived from the trained model. The Pareto front depicts that the exergy 
efficiency is ascending by increasing the heat cost function. In view of the 
fact that by increasing the exergy efficiency the energy efficiency is 
descending which cause to increase the cost function. There is an ideal 
point in Pareto front that never be accessible. The optimum point is 
obtained to the closest of the ideal point. At this point the heat cost function 
is equal to 0.015387 ($/kWh) and the exergy efficiency is equal to 10.01 (%). 


Investigation of aperture length and width disclosed that increasing the 
aperture length and width cause to increasing the exergy and energy 
efficiency. The cause of following content is extracting analysis of a 


household PTSC from a comprehensive trained model. Nevertheless, it is 
logical to select the aperture dimensions according to available commercial 
dimensions and constructability, so dimensions are selected to 2.44 (m) and 
2.29 (m) for the length and width, respectively. The concluded continuous 
and categorical variables are presented in Table 17. In addition, Table 18 
depicts related key parameters from analysis. 


Figure 10 depicts the variation of temperatures of HTF and tubes along the 
aperture length. The inner and outer temperatures of tubes is very close 
together, so in order to obvious illustration, it is showed just one of them. 
The outside surface of the receiver has the highest temperature, and the 
outer surface of the glass envelope has the lowest temperature. 
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Figure 10: The variation of temperatures along the length 


Table 17: Concluded continuous variables from optimization of the trained model 


N Variabl Uni 


Values 
O. e t 


1 D2 m 0.014 


2 R,-R; 
3 i 

4 Ww 

5 ih 

6 OT, 

7 P. annulus 
8 Xue 


9 Xx absorber 


10 


annulus 


0.0074 
2.44 
2,29 

0.31214 
60.7 

0 

Water 


Coppe 
r 


Argon 


Table 18: Other related key parameter that outcomes from analysis 


oF Variable Unit Values 
1 Ti i 62.79 
2 Du w 2685 
3 Energy % 74,33 
efficiency 
4 Exergy % 10.01 
efficiency 
fs) Exergy Ww 
destruction 3026 
6 AP Pa 7313 
7 Ryo - 62365 
8 Velocity m/s 2.053 


9 OF $/kKWh 0.015387 


4.5. variation calculation and sensitivity analysis 


Figure 11 illustrates the variation of energy efficiency, exergy efficiency, 
heat cost, energy CO, emission cost and exergy CO, emission cost by 
absorber inner diameter. By increasing the absorber inner diameter from 
0.01 cm to 0.08 cm, the heat cost increases. The exergy efficiency and 
energy efficiency are descending except at very small diameters. Due to 
high pressure drop and related increasing pump work, both efficiencies are 
decreasing at tiny diameters. By increasing the diameter, the efficiencies 
are declining due to decreasing in Nusselt number and convection 
coefficient that causes to decrease heat gain and increase exergy 
destruction. 
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Figure 11: The variation of energy efficiency, exergy efficiency, heat cost, energy-based C O7; 


emission cost and exergy-based C O; emission cost by absorber inner diameter 


Figure 12 illustrates the variation of energy efficiency, exergy efficiency, 
heat cost, energy-based CO, emission cost and exergy-based CO, emission 
cost by distance between absorber tube and glass envelope. By increasing 
the distance form 0.5 cm to 5 cm, the efficiencies are descending due to 
declining convection heat transfer coefficient at annulus space. The energy- 
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based CO, emission cost and exergy-based CO, emission cost are shown as 
a function of prior parameters. 
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Figure 12: The variation of energy efficiency, exergy efficiency, heat cost, energy-based carbon 
dioxide emission cost and exergy-based carbon dioxide emission cost by distance between absorber 
tube and glass envelope 


Figure 13 depicts the variety of energy efficiency, exergy efficiency, heat 
cost, energy-based CO, emission cost and exergy-based CO, emission cost 
by mass flow rate. By ascending mass flow rate, heat gain and destruction 
terms are ascending and the variation of destruction is greater than heat 
gain term, so the exergy efficiency is decreasing. Ascending the mass flow 
rate causes to increase the heat gain due to reach greater heat transfer 
coefficient, so the heat cost function is ascending. At this analysis, the 
energy efficiency is influenced by two primary parameters, the pressure 
drop and heat transfer coefficient. By increasing the mass flow rate both of 
them are increasing, but at low mass flow rates the convection heat transfer 
coefficient is dominant. Nevertheless, the energy efficiency is first 
ascending and then is descending. The energy-based CO, emission cost and 
exergy-based CO, emission cost are shown as a function of prior presented 
parameters. 
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Figure 13: The variety of energy efficiency, exergy efficiency, heat cost, energy-based carbon dioxide 
emission cost and exergy-based carbon dioxide emission cost by mass flow rate 
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Figure 14: The variation of energy efficiency, exergy efficiency, heat cost, energy-based carbon 
dioxide emission cost and exergy-based carbon dioxide emission cost by input temperature of fluid 


Figure 14 illustrates the variation of energy efficiency, exergy efficiency, 
heat cost, energy-based CO, emission cost and exergy-based CO, emission 
cost by input temperature of fluid. The energy efficiency is descending by 
increasing the input temperature, since the temperature difference between 
fluid and surrounding will be descending and subsequently the heat cost 
will be increasing. One of the exergy destruction reasons is temperature 
gap between hot and cold sources. Due to decreasing the temperature 
difference, the exergy destruction is decreasing that causes to elevate the 
exergy efficiency. 


Figure 15 depicts the variation of energy efficiency, exergy efficiency, 
heat cost, energy-based CO, emission cost and exergy-based CO, emission 
cost by annulus pressure between absorber tube diameter and glass 
envelope. Energy and exergy efficiency curves are coincided by increasing 
the pressure since their slopes are equal together based on related y-axis 
scales. By increasing the pressure, the heat transfer convection will be 
more substantial as heat loss. Increasing heat loss term in the efficiency 
relations causes to decrease them. It is obvious that a fluid with negligible 
pressure tends to eliminate heat transfer convection. 
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Figure 15: The variation of energy efficiency, exergy efficiency, heat cost, energy-based carbon 
dioxide emission cost and exergy-based carbon dioxide emission cost by annulus pressure between 
receiver tube diameter and glass pip 


Figure 16 illustrates the variation of energy efficiency, exergy efficiency, 
heat cost, energy-based CO, emission cost and exergy-based CO, emission 
cost by the length of collector. The efficiencies are increasing with the 
length. By increasing the length, ascending in heat gain is more than gained 
solar radiation. The heat cost is decreasing. The reason is that regarding to 
heat cost equation, by ascending the length of PTSC, the heat gain rate is 
increased more than the numerator of heat cost equation. The exergy 
efficiency is ascending, although the difference in temperature between 
heat sources and pressure term in exergy efficiency cause dissipating 
exergy, why the heat gain rate term is dominant over the mentioned terms. 
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Figure 16: The variation of energy efficiency, exergy efficiency, heat cost, energy-based carbon 
dioxide emission cost and exergy-based carbon dioxide emission cost by the length of collector 


Figure 17 shows the variation of energy efficiency, exergy efficiency, 
heat cost, energy-based CO, emission cost and exergy-based CO, emission 
cost by the width of collector. Similar to the variation of length, by 
increasing the length, the efficiencies are ascending and the heat cost 
function is descending. Furthermore, the reasons are the same based on 
thermal analysis in assessed bounds. 
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Figure 17: The variation of energy efficiency, exergy efficiency, heat cost, energy-based carbon 
dioxide emission cost and exergy-based carbon dioxide emission cost by the width of collector 


The energy efficiency, exergy efficiency, heat cost, energy-based CO, 
emission cost and exergy-based CO, emission cost are calculated and depicted 


to the various presented categorical variables in Figure 18, Figure 19, Figure 20, 
Figure 21, Figure 22. 
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Figure 19: The exergy efficiency for the various categorical variables 
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Figure 20: The heat cost for the various categorical variables 
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Figure 21: The energy-based carbon dioxide emission cost for the various categorical variables 
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Figure 22: The exergy-based carbon dioxide emission cost for the various categorical variables 


4.6. sensitivity analysis 


As a result of utilizing the proposed machine learning approach, a 
sensitivity analysis of all continuous and categorical variables is extracted. 
Figure 23, Figure 24, and Figure 25 depict the importance of each input on 
a percentage scale. According to sensitivity analysis of energy efficiency, 
exergy efficiency, and heat cost function, the most important variables are 
aperture width, input temperature, and inner diameter of the absorber. 
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Figure 23: Percentage of sensitivity for energy efficiency for each input variable 
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Figure 24: Percentage of sensitivity for exergy efficiency for each input variable 
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Figure 25: Percentage of sensitivity for heat cost function for each input variable 
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4. Conclusion 


The present study provided a prospective comprehensive approach to dealing 
with complex problems and obtaining quick responses from high-time solution 
problems by applying an appropriate machine learning method for academic 
research and industrial applications. The article includes an investigation of a 4E 
analysis and optimization of a complex renewable energy system with machine 
learning, which needs to solve multiple unknowns and multiple equations with non- 


linear terms. Through the proposed machine learning approach, assessing various 
inputs and variables to gain optimum conditions is provided. In this respect, after 
generating a comprehensive code for 4E analysis of a renewable energy system, 
validating with reference [52], obtaining proportionate inputs and outputs data and 
scaling, and providing a developed code to model an accurate machine learning 
model through the presented prospective approach, the considerable following 
consequences resulted: 


¥ A comprehensive and prospective detailed approach of machine 
learning associated with complex and high-time solution engineering 
problems was proposed and discussed in detail. 

¥ The proposed supervised machine learning approach utilized through 
a complex renewable energy system for a PTSC 

v After training the machine learning model from the proper inputs and 
outputs data under governing the presented non-linear constraints, 
dealing with analysis, obtaining quick solutions and feedback, and 
optimization with low time was gained and the problem transposed 
from a cumbersome complex problem to a simplified and convenient 
problem 

Y Utilizing the presented approach and the trained machine learning 
model for optimization revealed the best options for categorical inputs 
and continuous variables up to 1143 times faster process through a 
substantially more convenient and _ predictable method. The 
categorical variables including water for HTF, copper for absorber 
type and argon for annulus gas type and the obtained values for 
continuous inputs are 0.014 m for inner diameter, 0.0074 m for 
distance between inner and outer pipes, 0.31214 kg/s for mass flow 
rate, 60.7 °C for input HTF temperature and O pa for annulus gas 
pressure. The related optimum objective function resulted in 74.27 % 
for energy efficiency, 10.003 % for exergy efficiency, 0.0154 ($/kWh) 
7.07057 x10 *$/kWday for energy-based CO, emission cost and 
7.07048 x10* $/kWday for exergy-based CO, emission cost. 

¥ Furthermore, it is possible in a substantially more convenient and 
quick way for the simplified problem by machine learning model to 
perform post-processing. Thus, a comprehensive parametric 
investigations was illustrated for all input variables and objective 
functions in multi-scale depictions which assess the trends and 
turning points. 

Y In addition, a sensitivity analysis of variables is extracted from the 
trained machine learning model for objective functions, including 
energy efficiency, exergy efficiency and heat cost function. 


Nomenclature 


A Surface area of greenhouse (m7) | Abs Absorber 
D Diameter of pip Ch Chemical 
C; unit cost of fuel ($ kg") DE Direct energy 
IM je; | the fuel mass flow rate (kg s') | IDE Indirect energy 
C Unit cost ($/kg) CC Combustion chamber 
C, Specific heat (kJ/kg K) WF working fluid 
C Cost rate ($/s) Cond Condenser 
Cony amas’ at rate ($ h ‘ Exit 
Coo, Unit cost of casa dioxide ($ ae Ganersios 
f friction factor Int Intercooler 
activation function and the eae enue eonlecHon 
F hidden layers ourpute LW,,| weights ee the hidden 
ayer 
a, | feed the output layer's neuron | },,) gee ai e i 
Meco, aaa eee ere n Required sample size 
a | Distance ofebiectvencton | pr | Prandtel number 
eX Specific exergy (kJ/kg) Ex Exergy rate (kW) 
g Constant gravity Ex, | Exergy destruction rate (kW) 
f;, | Value of ith objective function 
Fa Partition function Subscripts 
h Specific enthalpy (kJ/kg) 0 Ambient condition 
| Interest rate In Inlet 
L Greenhouse length Com Compressor 
in Mass flow rate (kg/s) Out Outlet 
n System life time ph physical 
N Operating years Ref Reference 


P Pressure (kPa, MPa) SS Strong solution 
Q Heat transfer rate (kW) Tur Turbine 
required carbon dioxide rate 
R for a greenhouse per unit of WS Weak solution 
area 
Ss Specific entropy (kJ/kg K) Greek Symbols 
Nu Nusselt number 4 E Heat exchanger effectiveness 
k Conductivity Q Maintenance factor 
T Temperature (°C) Dax Equivalence ratio of fuel 
WwW Weight of units n Efficiency (%) 
Wi Greenhouse width p Energy function 
W Work production rate (kW) Abbreviations 
xX Mass fraction ORC Organic Rankin Cycle 
WA Cost ($) CRF Capital recovery factor 
Z Cost rate ($/s) HHV High Heating Value 
AH Change in enthalpy (kJ) GA Genetic algorithm 
AS Change in entropy (kJ/K) HE Heat exchanger 
HTF Heat Transfer Fluid 
CCH Combined Cooling, Heating, 
and Power 
GI Cost index 
re Deep Belief Neural Network 
GH Greenhouse 
GT Gas turbine 
LHV Low Heating Value 
AT temperature difference (°C) ARC Absorption Refrigerant Cycle 
COP Coefficient of Performance 
CRE Capital recovery factor 
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